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ABSTRACT 


An  important  consequence  of  the  quasi-linear 
theory  of  plasmas  is  that  a  weakly  turbulent  plasma  may 
be  regarded  as  a  gas  of  quasi-particles  such  as  plasmons 
and  photons  undergoing  collective  oscillations  and 
interacting  with  particles.  In  fact  this  process  may 
be  thought  of  as  an  emission  and  absorption  of  plasmons 
by  particles.  The  wave-wave  interaction  involves  only 
interactions  of  quasi-particle  among  themselves.  Some 
non-linear  plasma  phenomena  have  been  studied  by  starting 
with  the  Hamiltonian  and  making  use  of  the  second 
quantization  formalism. 

In  the  second  chapter  a  set  of  new  quasi-linear 
kinetic  equations  have  been  derived.  The  third  chapter 
contains  the  study  of  wave-wave  interaction  and  the  last 
chapter  deals  with  the  effect  of  wave-wave  interaction 
on  the  quasi-linear  theory. 
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A  BRIEF  REVIEW  OF  THE  NON-L INEAR  PLASMA  THEORIES 

INTRODUCTION 

Preliminary  studies  in  the  theory  of  non-linear  plasma 
have  been  made  by  Drummond  and  Pines  [l]  and  by  Vedenov, 

Velikov,  and  Sagdeev  [2  j .  Their  work  has  been  described  as 
a  quasi-linear  theory.  An  important  deduction  from  this 
work  is  that  certain  unstable  modes  persist  for  a  considerable 
length  of  time.  In  the  theory  the  particle-distribution 
function  is  split  into  two  parts  -  a  slowly  varying  part 
called  the  background,  and  a  rapidly  varying  part  which 
averages  out  to  zero.  This  rapidly  varying  part  represents 
a  system  of  oscillations  with  randomly  distributed  phases, 
whereas  the  slowly  varying  part  when  combined  with  the  linear 
terms  in  the  collisionless  Boltzmann  equation  yields  what  is 
called  a  quasi-linear  equation.  However,  when  the  non-linear 
terms  are  combined  with  the  rapidly  varying  part,  there 
results  coupling  between  various  modes.  By  averaging  over 
random  phases,  the  equation  for  the  slowly  varying  distribution 
is  obtained  in  the  form  of  a  diffusion  equation.  This 
diffusion  establishes  an  equilibrium  while  mode-mode  coupling 
causes  damping  of  the  spectrum. 

A  further  deduction  from  the  quasi-linear  theory  [1,3] 
is  that  a  weakly  turbulent  plasma  may  be  regarded  as  a  gas 
of  quasi-particles  such  as  plasmons  and  photons  undergoing 
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collective  oscillations  and  interacting  with  the  particles. 

In  fact  this  process  may  be  thought  of  as  an  emission  and 
absorption  of  plasmons  by  particles.  The  mode-mode  inter¬ 
action  involves  only  interactions  of  quasi-particles  among 
themselves.  Working  in  the  Hamiltonian  formalism  instead 
of  the  Boltzmann  equation,  permits  the  use  of  powerful  field 
theoretic  methods  which  yield  the  additional  spontaneous 
emission  term  in  the  kinetic  equations. 

Classically  one  has  to  solve  coupled  equations  for  one 
and  two  particle  distribution  functions  in  order  to  account 
for  this  very  term  which  is  essential  for  an  approach  to 
equilibrium  [4] . 

Since  the  quasi-linear  theory  deals  only  with  the 
interaction  between  particles  and  waves,  there  exists  a 
possibility  for  including  in  the  theory  the  interaction 
among  the  waves.  Some  estimates  of  these  effects  have  already 
been  made  [l,6]  for  the  potential  waves,  but  these  are  based 
on  the  Boltzmann  equation  formalism.  Asymptotic  perturbation 
theories  as  applied  to  hydrodynamic  equations  have  also  been 
used  [7,8] . 

The  study  of  non-linear  plasma  theories  is  necessary 
for  the  understanding  of  plasma  instabilities,  since  oscilla¬ 
tions  which  are  unstable  in  the  linear  approximation,  may  be 
stable  in  the  non-linear  case.  Langmuir  waves  are  most 
effectively  excited  in  a  plasma  under  certain  conditions  by 
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an  impinging  beam  of  charged  particles.  The  non-linear 
interactions  convert  the  Langmuir  into  transverse  oscillations, 
which  may  be  experimentally  observed.  As  a  rule,  the  generated 
oscillations  have  phase  velocities  of  the  order  of  the  beam 
velocity.  If  the  oscillations  can  be  diverted  by  non-linear 
effects  from  that  region  of  the  phase  space  in  which  they 
can  be  generated  by  the  beam,  with  sufficiently  high  efficiency 
then  cascade-like  development  of  the  stimulated  emission 
process  ceases  and  the  instability  is  stabilized.  This  example 
further  illustrates  the  general  premise  that  only  the  investi¬ 
gation  of  non-linear  effects  can  determine  the  practical 
significance  of  various  plasma  instabilities. 

The  non-linear  wave  interaction  that  leads  to  the 
Langmuir  oscillations  from  other  types  of  waves  is  therefore 
of  interest  for  plasma  confinement  problems. 

Another  example  of  a  non-linear  interaction  occurs  in 

the  earth’s  radiation  belts.  In  the  outer  radiation  belt 

4  6 

electrons  with  energies  of  5  x  10  -  2  x  10  ev  have  a 

— ■  2  —  f)  —  3 

measured  density  between  10  and  3  x  10  cm  .  This  density 

is  far  higher  than  can  be  expected  in  a  Maxwellian  distribution 

of  a  cold  plasma  in  which  the  average  particle  energy  is  1  ev 

3  3 

and  the  number  density  is  10  /cm  .  This  phenomenon  can  be 
explained  if  one  accepts  the  consequences  of  the  quasi-linear 
theory  which  predicts  the  possibility  of  a  statistical  accelera 
tion  of  particles  by  a  wave. 
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In  the  succeeding  chapters  quantum  mechanical  theory 
has  been  applied  to  investigate  wave-particle  and  wave-wave 
interactions.  The  quasi-linear  equations  thus  obtained  agree 
with  the  previously  known  work[l,2j.  The  plasma  model  con¬ 
sidered  neglects  the  collisions  between  particles,  and  no 
external  magnetic  field  is  supposed  to  be  present.  The  quasi- 
linear  equations  have  been  solved  analytically  and  exact 
numerical  results  for  the  steady  state  have  been  obtained. 
These  agree  well  with  the  results  obtained  numerically  by 
Drummond  and  Pines  [l | .  In  the  last  chapter  the  effect  of 
transverse  modes  on  quasi-linear  theory  has  been  studied. 


CHAPTER 
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QU AS I- L INEAR  PLASMAS  THEORY 


A  plasma  can  be  described  as  a  spectrum  of  collective 
oscillations  or  plasmons.  The  frequency  and  the  phase 
velocity  of  propagation  of  the  plasma  waves  are  determined  by 
the  dispersion  equation  which  contains  the  macroscopic 
parameters  of  the  plasma  such  as  temperature,  density, 
external  magnetic  field,  etc.  This  is  the  basis  of  the 
concept  that  all  the  particles  comprising  the  plasma  partake 
in  the  plasma  waves.  In  such  a  situation  the  dynamical 
aspect  of  the  particles  can  be  described  by  fluid-dynamical 
equations.  However,  this  being  a  collective  behaviour 
description,  the  phase  velocity  is  appreciably  larger  than 
the  individual  particle  velocities.  In  fact  one  is  implicitly 
assuming  so,  when  using  the  fluid  idynamical  equations. 

When  one  examines  the  damping  or  growth  of  waves  which 
is  determined  by  the  "fine  structure"  of  the  particle- 
distribution  in  phase  space,  the  role  played  by  the  resonance 
particles  becomes  important.  The  resonant  particles  are 
described  by 

_ y. 

a>.->  -  k . v  n  lOyj  n  0,1,2,  ...  (2.1) 

k  H 

where  cog-  and  £  are  the  frequency  and  the  wave  vectors 

characterizing  the  wave,  v  is  the  particle  velocity,  and 

eH 

a)  "the  cyclotron  frequency"  =  —  .  (H  =  external  magnetic 
H  rnc 

field,  m  =  mass  of  particle) .  Such  particles  are  capable  of 


- 
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exchanging  energy  with  the  wave  and  can  thus  damp  or  excite 
it.  We  shall  only  be  concerned  with  the  situation  when  n  =  o  or 
H  =  o  (no  external  magnetic  field) .  The  damping  of  the 
plasma  in  thermodynamic  equilibrium  is  found  to  be  dependent 
on  first  derivative  of  the  particle  distribution  function 
f  (v)  ,  v  =  iyk. 

A  linear  theory  is  unsatisfactory  since  the  growth  in 
amplitude  of  a  wave  is  unbounded.  A  quasi-linear  theory 
such  as  that  of  Drummond  and  Pines  [l] ,  predicts  a  process 
such  as  that  in  figure  (1)  in  which  a  particle  and  wave 
interact  giving  riseto  a°ne  particle  final  state. 


Fig.  (1) 


Fig.  (2) 


Fig.  (3) 
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1.  Particle-wave  Formation 


For  an  electron  plasma  with  zero 
field,  the  Hamiltonian  PI  is  given  by 


external  magnetic 


1  I 

2  i  <  j 


?  (P,-  ~ 

l 


2 

e 


e 

c 


it  (xi))2  + 


1_ 
8  TT 


(E2  +  H2)  dx  + 


(2.2) 


where 

m 


-> 


H 


-> 
x . 

l 

-e 

and  A 


mass  of  the  electron 
canonical  momentum 

electrical  field  for  the  transverse  waves 
magnetic  field  for  the  transverse  waves 
coordinate  of  the  particle 
charge  of  the  electron 
vector  potential 


The  summation  is  over  the  particles  and  the  gauge  chosen 
is  the  one  in  which  the  scalar  potential  is  zero.  The 
vector  potential  A  can  be  split  into  two  parts 

A  =  Al  +  At  (2.3) 

such  that  At  is  derivable  from  the  gradient  of  a  scalar 

• 

function,  and  the  divergence  of  AT  vanishes.  Thus  A^  can 
be  regarded  as  the  vector  potential  due  to  longitudinal 
waves  and  A^  to  the  transverse  waves.  However,  since  in 

this  chapter  we  shall  only  be  concerned  with  quasi-linear 


, 
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theory,  i.e.  particles  interacting  with  the  longitudinal 
wave,  we  shall  neglect  A^  for  now,  and  deal  with  a  non 
vanishing  in  the  work  of  the  next  chapter.  Let  us  make 
the  Fourier  expansion  of  A: 

A  (x ,  t )  =  E  A  (kV,t)  elK*X  (2.4) 

k 

since  A(x,t)  is  real: 


A  (■ 


•  )  = 


->* 

A 


(£,  t) 


(2.4a) 


For  k  <  k^  and  k  >  kD ,  the  Hamiltonian  is  separable  into 
two  parts  as  shown  by  Bohm  [  31  ],  where  kD  is  the  Debye 

wave  vector. 


II 


1_ 

2m 


2 

Pi 


+  V 


E 

k 


3^.  (k ,  t) 


8  TT  C 


3 1 


k 


D 


-*■  * 

3 A  (k  ,  t) 
3 1 


cop 


8ttc‘ 


A (k, t) 


*  A*(k,t) 


2mc 


E 

i  ,k 
k<kD 


•  A  -* 

-*■  ^  /i  j.  \  ik .  x  . 

.  .  A  (k ,  t )  e  i  + 

l  ' 


15 


+  'A  (k ,  t)  .  p. 


iit .  x . 

e  lj 


+ 


1 

2 


E 

i<  j 

k>  k 


x .  -x  . 
i  I) 


(2.5) 


D 


where  V  is  the  volurae,  it  stands  for  A.  .  The  second  term  is 

JLj 

the  Hamiltonian  of  a  harmonic  oscillator  with  coordinates  and 

V  3  z\* 

momenta  A(k,t)  and  — -  .  The  frequency  of  such 

A  ^  O  t 

4ttc 
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oscillators  is  evidently  co  .  The  third  term  is  due  to 

P 

the  interaction  between  the  particle  and  the  wave,  and  the 
fourth  is  due  to  particle-particle  effects.  For  such 
oscillators  we  can  define  according  to  Davydov  [33  ]  the 
normal  coordinates  given  by 


■/ 


y  wp 

2h 


AX^t )  + 


J  2yjW 


“>■  .  r>- 

n  (k,t) 


(2.6) 


r  u  ^p" 

/  2h 


* 

•~y  .  . 

A(k,t)  - 


/2uft 


0) 


n  (£,t) 


(2.7) 


where  y  = 


,  and  n  (£,t)  is  the  field  momenta 


4ttc‘ 


defined  above.  It  is  easily  shown  that  the  Poisson’s 
bracket 


{a  (k)  ,  a  (k '  )  } 


s  si?ic' 


(2.8) 


and  the  equations  of  motion  follow  from 


3a  (k) 
3 1 


=  {a  (k) ,  H}  = 


3  H 


3a  (k) 


(2.9) 


The  transition  to  quantum  mechanics  is  now  made  as  follows 


{a  (k)  ,  a*  (k ' )  }  =>  -  i  (a(k),  a*(k')>  =>  -  1 


'h 


a  (k) 


a*  (k ' )  -  a* (k')  a  (k) 


; 


=  -  y-  fiyy  i 
'h  kk ' 


(2.10) 
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so  that,  the  commutation  relation  is 

[a<k>,  a*(k')]  = 


(2.11) 


+ 

Thus  in  the  transition  to  quantum  mechanics  a(k)  and  a  (k) 
become  the  second  quantized  operators  for  plasmons  of 
momentum  hk.  If  |  N  (it)  >  is  a  state  vector  with  N  (it) 
plasmons  in  the  state  it,  it  follows 


+ 


a  (k)  j  ....  N (k)  . . . . >  = 


/N (k) +1 


....  N  (k ) +1  ....> 

(2.12) 


a  (it) 


n  (it) 


_ N(it)-i 


.  > 


(2.13) 


+ 

The  operators  a  (k)  and  a  (k)  create  and  destroy  respectively 
one  plasmon  from  a  state  k. 

The  second  quantization  Hamiltonian  for  the  wave- 
particle  interaction  can  simply  be  written  by  using  (2.6) 
and  (2.7)  and 

dx  l+(x,t.)  £  [  p..it(k,t)  + 

i,k  1 

^(}t,t)  .pj  elv’'^  l(x,t)  (2.14) 


II 


co-p 


e 

2mc 


where 


I  (x,t) 

=  £ 

b  (k) 

K  (x) 

(2.15) 

k 

+  ->  . 

I  (x,t) 

=  £ 

b+(k) 

+  -> 

V  (x) 

(2.16) 

k 
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where  £  (x)  is  the  plane  wave  function  for  the  particle, 

b+  and  b  are  the  creation  and  annihilation  operators 
for  the  particles. 

After  using  (2.6)  to  (2.16),  the  Hamiltonian  finally 
becomes  (after  neglecting  particle-particle  interactions.) 


H 


H  +  H .  . 

o  mt 


(2.17) 


with 


H 


i  a)  (a  a+  +  a+  a)  (2.18) 

z  k  p 


H.  , 
mt 


=  V  ’  eh  /  /  >h 
/  ,  1  me  W  2yu 

K^' 


P 


a  (k)  b+  (it1) 


b(ic'-S)  +  a+(ic)  b+(ic'-ic)  b(it') 


+  H.C 
(2.19) 


will  give  us  the  probabilities  for  emission  and 


absorption  of  plasmons  by  the  particles. 


2.  The  Development  of  the  Distribution 


If  a  larger  number  of  plasmons  are  emitted  in  a  state 
than  are  absorbed  from  it,  the  wave  is  unstable.  The  wave 
is  stable  if  the  reverse  happens.  We  have  to  calculate  the 
net  number  of  plasmons  in  a  state.  If  N  (ic )  and  f  (it)  are  the 
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number  of  plasmons  and  particles  in  a  state  k ,  then, 
schematically 


dN  (}t) 
dt 


or 


E 

it' 


^  _TV-£-  It' 


V» 


dN  (it) 
dt 


6  (hu>  +  Ek  , 


E 


k+k 


(2.20a) 


where  M  is  the  vertex  part  and  E's  are  the  particle 
energies . 

Similarly 


df  (it'  ) 
dt 


E  | M  |  2  N  (it) 
li  k 


h 

m 


2 


m 

h 


)2k2 


3 

3k  ' 


hit,  it 1 

m 


x 


3  f  (k  1  ) 

IF* 


+ 


E 

it 


3 

3k' 


it.  hit1 

m 


(2.21) 
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The  second  term  above  is  the  spontaneous  emission  term, 
which  is  absent  in  the  classical  theories. 

i  \2 

If  we  take  the  expectation  value  of  |  A^  |  between 
states  that  are  eigen  functions  of  the  number  operator, 
then 

|A(k)  I  2  : 


(k) 


y  w 


Therefore 


c 


to  2 

JP 

8n 


A  (k ) 


^p^  4iN  (k ) 

8  n  y  oo 


(2.22) 

where  's’  is  the  electric  energy  density  in  the  waves. 

The  classical  limit  of  (2.20)  and  (2.20a)  give,  the 
fundamental  equations  of  quasi-l.inear  theory.  If  one 
neglects  the  small  spontaneous  emission  term,  then  these 
equations  are 


i£  - 
at 


a  (v) 


3f  (v) 
3v 


(2.23) 


a  f  (v) 
at 


av 


a  ,  o  /  \  -t?  /  i  \  a  f  (v)  ) 

(  3  (v)  (5  (v,t)  — ^ — ) 


v;here  a  (v) 


.222 

4  t\  e  v 

moo 


and 


3  (v)  = 


3v 


8TfLe‘ 

2 

m  v 


(2.24) 


where  L  is  the  side  of  the  cube  containing  the  plasma. 
These  equations  are  identical  to  those  obtained  in  [l,2]. 
By  some  algebraic  manipulation  these  equations  yield 


3_ 

at 


[ 


-  1 


t 


3]f 

a  t 


a  (v) 


32 

av2 


(ft)]  =  0 


(2.24a) 
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Therefore 


o 

G 


at 


a  (v) 


3  v 


2-i)  =  c 

a 


(  -- C  )  =  constant 


(2.25) 


at  t  =  o,  f  =  F,  j  where  F-,  is  given  by  [l] 


Fo  = 


n 


Exp  [  -h  (  ~  )2  ] 


(  2  TT  )  hv 


+  4  x  10 


-4 


Exp 


[(  £  -  5)"  ] 


(2.26) 


and  c  =  C'  o  then 


c 


3  € 
3 1 


a  (v) 


3  2 
v 


a  ~P 

< \C  ) 


=  a  (v) 


3 

3v 


''v'l 

Fo  ^v)  - 


3  v 


(  iH 


a 


o  ) 


(2.27) 


In  the  steady  state  If-  is  described  by  the  equation 


"jS 

C  -  h, 


a 


Fo  (vj  dv 


v 


o 


Fo (vj  dv 


(2.28) 


or 


o  1 


e 


<f  ,3 


/■,- 


o 


Fo  (v)  dv 


v 

o 


Fo (v)  dv 


(2.28a) 


where 


Ly 


#-v  -p  (l) 

2L  p 
Z2 

v  Trm 


(t-  <;o) 


The  integrals  of  the  type 


v 

o 


Fo(vj  dv  can  only  be 


/■  i 

solved  numerically  [19,  15 J .  The  fig.  (4  )  shows  g 
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against  (  —  ) .  This  determines  the  steady  state  spectrum 
v 

of  electrostatic  energy.  It  is  worth  noting  that  when 


it 

3t 


o,  then 


9f 
3  v 


=  o. 


Thus  the  distribution  becomes 


flat  near  v 


w 


at  very  large  time  wpt  -> 


/k 


The  curve 


in  fig.  (4)  is  in  fairly  good  agreement  with  the  previous 
work  [l,  2,  3]  for  values  of  the  order  of  t  ^  sec. 


The  spectrum  C/  (  — —■  )  that  has  been  obtained  by  the 

methods  of  this  chapter  produce  a  particle  distribution 

which  differs  from  the  Maxwellian  distribution  only  for 

v  ^  wp  .  Details  of  the  other  interactions  mentioned  are 
k 

contained  in  the  Appendix. 
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COUPLING  0 F  THE  LONGITUDINAL  TO  TRANSVERSE  WAVES 

While  early  studies  of  non-linear  wave-wave  inter¬ 
actions  were  carried  out  by  Ginzburg  [  1 6 |  and  Bailey  [ 17  j 
it  is  only  more  recently  that  the  theory  has  been 
developed  to  such  an  extent  that  the  physical  mechanisms 
are  now  well  understood.  The  rapidly  growing  experimental 
research  in  this  area  of  research  has  provided  an  added 
incentive  for  the  pursuit  of  theoretical  studies.  To 
date  attention  has  been  directed  mainly  to  cases  of  weak 
non-linearity  wherein  the  wave-wave  modes  which  result  from 
the  weak  non-linear  interaction  can  develop  to  a  significant 
wave-amplitude  as  the  system  approaches  equilibrium. 

The  first  study  of  non-linear  interaction  specifically 
with  respect  to  Langmuir  oscillations  was  due  to  Sturrock  [~18  | 
Further  development  of  the  theory  was  made  by  Drummond  and 
Pines  [l] ,  and  Zheleznyakov  [19] .  The  first  attempt  to 
study  the  non-linear  effects  in  the  framework  of  kinetic 
approach  was  made  by  Kodomtsev  and  Petviashvili  [20 j, 

Karpman  [21 | ,  Galeev,  Karpman,  and  Sagdeev  [22],  Camoc  [23] 
and  Rudakov,  Vedenov  et  al  [24] .  Some  non-linear  problems 
connected  with  Langmuir  waves  were  also  studied  by  Silin  [25] 
who  used  the  correlation  method  developed  by  Bogolyubov  [34], 
This  approach  is  mathematically  so  involved  that  it  does 
not  permit  extensive  use  in  some  difficult  problems. 


" 
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In  this  chapter  only  the  transverse-longitudinal 
wave-wave  interaction  will  be  considered  in  detail.  The 
problem  of  excitation  of  transverse  plasma  modes  [27,29] 
and  the  energy  transferred  to  them  [26]  by  non-linear 
interaction  of  longitudinal  waves  has  previously  been 
studied . 


1 .  Wave- Interaction  Formalism 

For  an  electron  plasma  without  any  external  field, 
the  Hamiltonian  H  is  the  same  as  given  by  (2.2)  together 
with  the  definition  (2.3) 


1  I 

2  i  <  j 


A  (x±)2)  + 


1_ 

8tt 


(E2 


+ 


H2) 


dx  + 


(3.1) 


The  Hamiltonian  (3.1)  can  be  separated  into  two  parts,  the 
collective  and  the  particle. 


H  = 


I 

i 


2 

e 

1  2 
2mc 


+ 


1_ 

8  TT 


(E2  +  H2)  dx  + 


1_ 

8  TT 


(3.2) 


V7here  E  is  the 
J-J 


longitudinal  electric 


field . 
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-V  -V 

In  (3.2)  the  (A.p)  term  which  is  responsible  for 
particle-wave  interaction  is  absent  unlike  the  case  of 
second  chapter.  This  is  so  because  for  large  wavelength 
plasma  oscillations  only  a  few  particles  are  resonantly 
coupled  with  the  wave .  However,  in  the  Hamiltonian  the 
contribution  due  to  transverse  waves,  which  was  neglected 
in  the  formulation  of  Chapter  2,  is  now  a  significant 
factor.  Expanding  the  vector  potential  as  before  in  a 
cubic  box  of  length  1 

t  r  4-\  £  t  /,  tv  ik.x 

^L,T(Xrt)  =  k  AL,T(t'k)  e 

with 

A*L/T(j"'t)  =  AL,T("j"'t)  (3‘4) 


Using  the  Poisson's  equation,  the  Hamiltonian  becomes 


H  =  H  +  H  +  H.  ,  , 

T  L  int  ' 


(3.5) 


where 


->  • 


HT  =  V  1  t  2  3 1 

k  Sue 


1  3A>j<p  x  3  A  „  .  k 


(k,t).  ^~T  (k,t)  +  At  (k,t) . 


3  r 


-> 


at  (it ,  t)  + 


— ^2  ?T(S,t)  .A*  (ic,t))  , 
8ttc 


(3.6) 


H.  =  V  s  |_L  (k,t)  .  (k,t)  +  -A 


r>  0  2  3 1 

k  8ttc 


8ttC 


2  L 


_y*  -> 


,Al  (k,t) ) , 


(3.7) 
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H.  . 
lilt 


WP  e  y 

o  3  V  U 
m8  tt  c  kk 


XL(jc,t)  .AL(kv'  ,t)  +  AT(lt!  ,t)  .AL(]t,t)  + 


-> 

A, 


,  ~y 


(k,t) 


A. 


(k  '  ,  t )  + 


at  (£' 


/t) 


.AT(kft) 


A  (k-!  k  1  ,  t ) 


— ^ 

/  V  i  ] 
\Atj 


it  i  ^ 


and  V  ,  the  volume. 


(3.8) 


H,^  is  the  Hamiltonian  for  the  transverse  waves, 

Ht  for  the  longitudinal  waves  and  H.  ,  for  the  inter- 
action  between  the  two.  That  part  of  the  Hamiltonian 
which  is  responsible  for  the  self  interaction  of  the  waves 
has  been  neglected. 

For  the  purpose  of  forming  the  plasmon  Hamiltonian, 
it  is  convenient  to  introduce  the  canonical  variable  for 
the  longitudinal  and  transverse  waves, 


v 

uL,t  .  2  at  L,T 

'  4ttC 


(3.9) 


*l,t 


% 


L,T 


(3.10) 


where  tt  is  the  momentum  and  $"  the  coordinate.  Next  v/e 
can  form  the  plasmon  normal  coordinates  defined  by 


L,T 


'  2  -h  L  ,  T 


,2u-hwT  L  (k ) 


L , T ,  (3.11) 


L,  t 


/uw  T  (k)  * 

/  _L'it _ _  a 

2h  L ,  T 


'2yhiwT  L  (k) 


* 

TT 


L,T'  (3.12) 
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where 


y 


v 


4ttc‘ 


to., 


2 

L 


2  ,  2  n  2  2,  2, 
to  and  10^=  (k  c  +  to  ) 
P  T  p' 


(3.13) 


The  a's  and  a  ' s  are  the  annihilation  and  creation 
operators  of  second  quantization  for  the  plasmons 
Using  (3.9-13)  we  can  write  (3.6-8)  as 


H 


~  di  £  toT(k)  (aTa*-fa*aT) 


(3.14) 


H, 


1  ^  E 

2  *  k 


“L  (aLaL+aLaL> 


(3.15) 


H.  ,  =  H,  +  Hn 

mt  1  2 


(3.16) 


H 


to 


Pe 


ri 


3/2 


^ 

-  TV  2  / - o - 7TT  (k+k1)  cos0  [am  (k)  ar  (k  '  ) 

1  m8^c3  V  kk'  /Vr-8coT(k)y  L  T  L 


aL  (k+k ' )  +  aT  (k ) AL  (k 1 ) aL  (k+k ' ) +aT  (k ) aL  (k ' ) aL  (k+k 1 ) +aT (k  r ) 


aL  (k) aL  (k+k 1 )  +  aT  (k) AL  (k ' ) aL (k+k 1 )  +  H.C. 


(3.17) 


H. 


•to 


Pe 


•fi 


3/2 


V  E 


(k  *  )  com  (k)  0J„y 


K  cos9 


m8 tt c  kk  ’  /y  SooT 


(a*  (k)a*  (k ' ) aL  (k+k* ) +aT  (k)aT  (k1 ) aL  (k+k' ) +aT (k)aT (k ' ) aL (k+k’ ) 


+  aT  (k )  a^  (k  '  )  aL  (k+k  1  )  )  +  H.C, 


(3.18) 
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Hm  and  H  are  the  oscillator  energy  of  the  transverse  and 
longitudinal  waves.  is  due  to  two  longitudinal  waves 

colliding  producing  a  single  transverse  w ave  and 
represents  two  processes:  (a)  two  transverse  waves  collide 
to  produce  one  longitudinal  wave  (b)  one  longitudinal  and 
one  transverse  wave  interact  to  result  in  one  transverse 
wave . 

Thus  from  H.  we  can  write  the  matrix  elements  for 
mt 

the  respective  transitions.  Knowing  this  and  using  the 
"Golden  Rule"  for  time  dependent  transitions ,  we  will  be 
in  a  position  to  write  the  time  development  for  the  plasmon 
distribution.  * 

2 •  £ lasmon  Distribution 

If  N  (k)  and  n(k)  are  the  number  of  plasmons  for  the 
longitudinal  and  transverse  waves  respectively  and  if  the 
curly  line  denotes  longitudinal  wave  and  the  straight  line 
transverse  wave.,,  then  schematically 


' 
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Similarly  one  can  write  an  equation  for  (k) . 

% 

So  long  as  the  longitudinal  wave  vector  k  <<  —  the  last 

two  interactions  in  (3.19)  do  not  conserve  energy.  Thus 

in  this  limit  we  are  left  with  only  the  first  process. 

kL<<w  /c  implies  that  we  need  little  longer  wavelengths 

than  are  necessary  to  satisfy  X<<AD«  However  if  the  last 

two  processes  have  to  play  any  significant  role,  the 

particle-wave  interactions  cannot  be  neglected. 

In  the  limit  k  <<oj  /c 

L  p7 

dNOO  =2 *  JZ  I M.  |  2  [n(it’)  {  N  (Jc )  +  N(S'-it)  }  - 

dt  *i  k,  1'  L 

N  (it)  N(g'-ic)]  6  [h  ((0T(itI)  -  2w  ) ]  (3.20) 

=  y  |M'  |  2  fN(kV)  N (£'-£)  -  n(lt')  x 

dt  ^  ^  1 

{  N  (3c)  +  N  (£'-£)}]  6  fh  (ojT(it')  -  2u>  )]  (3.21) 


where  M-,  and  M~  are  the  matrix  elements  for  the  respective 
transitions  and  can  be  computed  simply  from  the  Hamiltonian 

i  i  2 

.(3.17).  If  we  take  the  expectation  value  of  | T| 
between  states  that  are  eigen  functions  of  the  number 
operator,  then 


< 


Vk) 


2 

> 


&N  (k) 


y  0) 


P 


(3.22) 
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<  AT(k)  > 


hn  (k)  _ 

pt  (k) 


(3.23) 


for  N  (k) ,  n  (k) >>  1 . 


We  shall  assume  that  the  longitudinal  spectrum  [32 J  is 
zero  outside  the  interval  AK  and  constant  within  AK  and 


2  2  tt  3 

of  magnitude  |E(t) |  ,  where  AK  =  — j  A^ 


If  at 

particular  instant  the  interaction  comes  to  life  and  if 
the  longitudinal  spectrum  is  known  at  that  time,  then  the 
transverse  energy  density  Wr  (t)  can  simply  be  written  down 
from  the  equations  (3.20)  and  (3.21). 


wT(t) 


32 


JL 

tt  a 


E 


o 


1  - 


a+  3~a  Exp [  - 


c  lr  I  2 

— 2  I Eo 1  T 
0) 

P 


3  ] 

(3.24) 


where  t(>tc)  is  the  interval  of  time  elapsed  after  the 

longitudinal  waves  had  started  interacting  to  generate 

transverse  waves,  t  being  the  characteristic  time  of 

c 

interaction . 


3  = 


8tt3/3  L6 


2tt 


4 

U) 


K 


m 


D 


a  = 


2,-42 
tt  V  3  co  e 
_ _ P _ 

2  7 
m  c 


L__ 
2  TT 


E 


in  the  absence  of 


i  2  v  T 

E  =  - — o  is  the  value  of 

°  2L3 

interaction,  x  being  the  Boltzmann  constant.  Therefore 


the  characteristic  time  for  the  growth  of  transverse 


■ 
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energy  is  given  by 


or 


T 


16/3tt 
6nl  2 


5 


(3.25) 


where  6n  is  the  change  in  the  electron  density  n  and  v 
the  thermal  velocity. 

If  one  calculates  the  rate  of  growth  of  energy  in 
the  transverse  modes  for  a>t<<l. 


dWT(t) 

dt 


3.2  x 


nxT 

(nXD)d 


(3.26) 


This  agrees  with  the  previous  work  [26] .  Thus  the 
present  formalism  makes  it  possible  to  derive  the  previously 
known  results  as  a  special  case.  Moreover  it  throws  a 
light  on  other  possible  third  order  interactions  and  how 
they  can  be  eliminated.  Another  contribution  of  the  present 
chapter  has  been  the  calculation  of  the  characteristic  time 
for  the  energy  transfer  to  transverse  waves. 
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CHAPTER  4 
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EFFECT  OF  INTERACTION  BETWEEN  TRANSVERSE  AND  LONGITUDINAL 

WAVES  ON  QUASI -LINEAR  THEORY 

The  result  is  developed  in  Chapter  2,  that  in  the 
quasi-linear  approximation,  the  plasma  oscillations  have  a 
steady  state  particle  distribution  function  with  a 
plateau  near  v  =  ~  .  To  date  no  work  has  been  done  on 
the  effect  of  non-linear  interactions  between  longitudinal 
and  transverse  waves  on  the  particle  velocity  distribution. 
However  it  is  important  to  examine  this  problem  since,  for 
times  which  are  sufficiently  long  for  the  steady  state  to 
be  attained,  wave-wave  interactions  become  important. 

By  similar  arguments  to  those  which  were  used  to 
derive  (2.20)  and  (3.19),  the  particle-wave  and  wave-wave 
effects  can  be  combined  into  one  single  equation.  As  before 
the  interaction  in  which  longitudinal  waves  decay  into ' two 
transverse  waves  will  not  be  considered.  Therefore 
schematically  one  obtains 
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The  first  summation  is  due  to  the  transverse- 
longitudinal  wave-wave  and  the  second  due  to  wave-particle 
interactions.  Upon  applying  the  principle  of  detailed 
balancing  for  time  independent  transitions  to  (4.1)  one 
obtains 


dN(k)  2tt 


Z  |M-,|Z  [N(jt)  {n  (ic '  )  -  N(k'-it)}  + 


->  a- , 


dt 


*  P 


N(it'-ic)  n  (it  ’ )]  6  [h(uT(ic')-2u  )] 
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Z 

& ' 


m2|2  [f(ic+ic') 


-  f(k')]  N  (k)  x 


4  ft  y  +  Ejj,  -  E^.'J 


(4.2) 


where  n(k)  is  number  of  transverse  plasmons  or  the  quasi¬ 
photons.  Similarly  one  can  write  for  n(it) 


dn  (k  1  ) 
dt 


2tt 

“ft 


Z 
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n(it') 
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(4.3) 


The  equation  for  df(k’)  is 

dt 


df  (k ' ) 
dt 
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+  I  lM2l2  k  Tk'  s  <“P  f  <k'>  > 


After  some  algebra,  the  fundamental  equations  for  the 
elementary  processes  become 
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at  av  (B(V)^  H> 


3  f  3 


(4.6) 


S.  and  ^ 


where  C2  and  are  the  electric  energy  densities  of  the 

L  J- 


longitudinal  and  transverse  and 
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v 
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a  (Y) 


3  (v) 


.222 
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(4.8) 


(4.9) 


In  the  steady  state  equation  (4 . 5^)  yields 
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A. 


(ck  +  /3  wp) 2  5  l  (k) 


dk 


’(ck  +  /3  “p)2['5L(k)  +~S-L( k  -  /3  jg.) 
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dk 


(4.10) 


If  k  >>  /3  jp 
c 


and,  of  course,  less  than  kD,  then 


(4.11) 


Substituting  (4.11)  into  (4.4),  we  see  that  the 
resulting  equation  together  with  (4.6)  are  the  equations 
of  the  quasi-linear  theory.  Thus  the  transverse-longitudinal 
wave  interaction  leads  to  a  redistribution  of  the  energy 
balance  in  the  various  modes,  the  spectral  shape  of  the 
longitudinal  modes  is  not  affected  by  this  type  of  interaction. 
The  distribution  function,  however,  displays  the  same 
characteristic  properties  of  Chapter  2. 
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APPENDIX. 


consider  the  term 


1_ 

2m 


<*i> 


(A.l) 


After  Fourier  analysing  (A.l)  becomes 
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A  (St)  A  (it') 
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This  further  is  equal  to 


H 

c 


1_ 

2m 


k, k 1 k" 


6  n  (it"  )  a  (it)  A  (it) 


i  (it+it'+it")  x.) 
e  i 


dx  (A. 3) 


where  6  n(k")  is  the  perturbation  in  the  electron 

density . 

This  is  non  zero  is  k"  =  -  (k+k1).  If  k"  ^  -  (k+k1), 

then  H  remains  the  same  as  in  (A. 2).  That  is,  it  is  not 
c 

possible  to  eliminate  the  particle  variables.  Therefore 
transforming  to  second  quantization  formalism  (A. 2)  becomes 
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+  H.C . 
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Thus  ore  can  see  that  H  gives  rise  to  the 

c 


interactions  mentioned  above. 
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